Multi-photon quantum interference in a multi-port integrated photonic device 
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Increasing the complexity of quantum photonic devices is essential for many optical information 
processing applications to reach a regime beyond what can be classically simulated, and integrated 
photonics has emerged as a leading platform for achieving this. Here, we demonstrate three-photon 
quantum operation of an integrated device containing three coupled interferometers, eight spatial 
modes and many classical and nonclassical interferences. This represents a critical advance over 
previous complexities and the first on-chip nonclassical interference with more than two photonic 
inputs. We introduce a new scheme to verify quantum behaviour, using classically characterised 
device elements and hierarchies of photon correlation functions. We accurately predict the device's 
quantum behaviour and show operation inconsistent with both classical and bi-separable quantum 
models. Such methods for verifying multiphoton quantum behaviour are vital for achieving increased 
circuit complexity. Our experiment paves the way for the next generation of integrated photonic 
quantum simulation and computing devices. 



Realizing quantum-enhanced information processors 
for tasks such as simulation and computation demands 
experimental systems of sufficient complexity that their 
dynamics cannot be efficiently determined using classical 
processors. Reaching this regime in practice, however, 
remains a critical open challenge. 

Integrated quantum optics provides great promise for 
enabling photonic experiments, which are otherwise gen- 
erally limited to relatively small-scale experiments, to 
reach this new regime of complexity. Chip-based fab- 
rication enables sophisticated networks involving multi- 
ple interfering pathways in a compact and stable phys- 
ical architecture, and pioneering work has shown the 
viability of this approach for the manipulation of the 
quantum properties of photons ^-[5]. To date, exper- 
iments have demonstrated up to three-photon higher- 
order terms from a single nonlinear photon source being 
coupled into the two input modes of a single interferome- 
ter [5], or alternatively, two-photon correlations in up to 
21 waveguide modes [HOU]. Building a photonic system 
capable of truly outperforming classical processors, how- 
ever, can only be achieved by simultaneously increasing 
the number of modes and interference nodes in the circuit 
and the number of photons distributed among them. 

There are two key outstanding challenges associated 
with this task of scaling up integrated photonic circuits 
to these larger systems. First, photon loss exponentially 
limits the complexity achievable in a quantum circuit, 
both in terms of the number of circuit elements and 
the number of photons that can be used effectively. In 
integrated photonics, significant losses arise from inter- 



facing the circuit with both sources and detectors and 
become more pronounced with increased photon num- 
bers [TT1 [T2] . Ultimately, losses are fundamentally lim- 
ited by the intrinsic optical properties of the medium 
and these clearly scale with the circuit size. Secondly, 
the monolithic nature of integrated architectures means 
it is also difficult to verify that the chip meets the design 
specification. In particular, it is not possible in general 
to access individual components in situ using the exter- 
nal input and output ports, nor is it always possible to 
configure ancillary access ports for injecting probe states 
or performing detection locally. On the other hand, ex- 
isting process tomography techniques for verifying the 
quantum operation of a full chip 1 1 31 114] become imprac- 
ticable once it becomes sufficiently complex. Instead, 
other simpler ways to measure the chip's transformation 
are required. 

In this work, we demonstrate a critical advance in the 
complexity of quantum integrated photonic devices by si- 
multaneously increasing the number of photons and the 
number of spatial modes used. Using a circuit in which 
three photons are distributed over eight modes and three 
coupled interferometers [I51U7| . we certify quantum op- 
eration beyond both the classical limit and what can be 
achieved with two photons using a hierarchy of higher- 
order photon correlation functions. As part of this, we 
also provide the first on-chip demonstration of a Hong- 
Ou-Mandel-type interference effect with three individual 
input photons. 

We further show that a critical step in verifying the 
correct quantum operation was characterizing the oper- 
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FIG. 1. Schematic of the multi-port waveguide circuit. 

The circuit consists of eight modes labelled a to h, ten beam 
splitters labelled rji to rjio, and three variable phase shifters 4>i 
to 03. Single photons are launched into a subset of b to g and 
the output of each mode was monitored by a single-photon 
detector. The ancillary modes a and h were not accessible for 
coupling. 

ating parameters of individual circuit components, and 
we introduce a simple loss-insensitive method to achieve 
this using classical light scattered from the device in the 
transverse direction. Together with knowledge of the cir- 
cuit topology, this method allows the full unitary trans- 
formation implemented by the device to be reconstructed 
and used to verify quantum operation. In this paper, 
we verify three-photon interactions which achieve a com- 
plexity sufficient to realise a next generation of on-chip 
quantum information protocols such as cluster-state gen- 
eration and teleportation. 

RESULTS 

The multiport waveguide circuit used in these exper- 
iments consists of three coupled Mach-Zehnder inter- 
ferometers (MZI) spanning eight spatial modes, with 
phase control inside the interferometers implemented by 
thermo-optic phase shifters (see Fig. [I] and Methods for 
a detailed description of the experimental apparatus). 
Our main aim is to show genuinely quantum operation 
of the circuit in a context which demonstrates its full 
complexity in terms of simultaneously increased number 
of independent input photons and number of interacting 
modes and interferometers. To do this, we inject individ- 
ual single-photon states into one mode of each of the in- 
terferometers (modes c, d and /) and measure the visibil- 
ity of three-photon nonclassical interference at different 
combinations of output ports. If the observed visibility 
is stronger than that predicted when the single-photon 
inputs are replaced with classical light, this acts as a 
witness to the desired quantum behaviour. 

To calculate the predicted classical bounds, we must 
assume that the circuit itself operates completely coher- 
ently (unitary operation) , so we first characterise the cir- 
cuit using an element-wise, loss-tolerant approach with 



classical input states, in the process avoiding any need 
for resource-intensive quantum process tomography. We 
confirm the reliability of our characterisation by first ob- 
serving two-photon quantum interference and comparing 
its behaviour with both a classical and a quantum model: 
the former is clearly inconsistent with the experimental 
results, while the latter shows good agreement. Finally, 
extending this method, we show that the observed three- 
photon interference measurements exclude with high con- 
fidence levels both the classical model, as well as quan- 
tum models involving biseparable states. The nature of 
the observed three-photon interferences, combined with 
the known topology of the circuit relative to the input 
photons, suggest that the results cannot be explained by 
using a simplified or restricted subsection of the illus- 
trated circuit. This implies that the experiment utilises 
the full available complexity of the circuit. 

Multiphoton integrated-optics experiments set strin- 
gent demands on performance with regard to photon 
loss pTJ [12]. Particular care needs to be taken to op- 
timise all experimental efficiencies, especially in experi- 
ments utilising down-conversion photon sources, such as 
this one, to minimise higher-order noise terms |18j . In 
order to demonstrate high-brightness multiphoton states 
"on-chip", we have combined a range of technical solu- 
tions for optimising the loss properties, including efficient 
pair-source heralding, optimal coupling of six-channel 
fiber arrays to the chip and use of a low-loss integrated 
platform (see Methods). These measures were critical to 
reaching the level of on-chip complexity achieved in these 
experiments. 

Characterising circuit operation 

The three coupled interferometers in our waveguide 
circuit, fabricated by UV direct-write technology on a 
silica-on-silicon substrate [SJ, involve ten beam splitters 
and three thermo-optic phase shifters, a circuit which 
has only previously been realised directly in a simpli- 
fied polarisation-based encoding with bulk optics [ToT - 
I17j . Temperature control by means of a thermo-electric 
Peltier element maintained stable beam-splitting ratios 
and phase offsets over many weeks. Individually char- 
acterising these parameters permits us to simulate the 
required multiphoton interference visibilities. 

We measured the beam splitter reflectivities, r]i to 7710, 
by sending continuous- wave laser light through each split- 
ter in turn and calculating the reflectivity using a ra- 
tiomctric analysis which is independent of coupling and 
transmission losses |12j . This technique uses four mea- 
surements for each beam splitter, coupling light in turn 
into each input port and recording the power at each 
output port. Due to the complex circuit topology, it is 
not generally possible to independently access the input 
and output modes, so instead, light scattered out of the 
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FIG. 2. Classical characterisation of integrated circuit 
elements, (a-b) CCD images of light scattered in the trans- 
verse direction from the waveguide. Red arrows denote the 
direction of propagation of the light. The white box denotes 
the integration region at the output for the ratiometric analy- 
sis, (c) The integrated intensity is used as part of a ratiome- 
teric analysis to determine r/s = 0.55 ± 0.02. (d) Interference 
fringe produced by scanning tf>i. The parameter values ob- 
tained from the theoretical fit include the zero-voltage phase 
offset. 



chip surface was used, to measure the output powers. The 
ratiometric calculation is insensitive to different scatter- 
ing efficiencies in the same way as to different interface 
coupling efficiencies. Figures [2ja) and (b) show a typi- 
cal example, with 100 mW of laser light coupled into the 
chip (note that splitters 4 and 7 only had one available 
input, since modes a and h were not accessible for cou- 
pling, see Fig. [T]). The input polarization was adjusted 
to maximize the amount of the transverse scatter, which 
was imaged using a CCD with a highly linear response. 
Integrating over a specified region then provides the re- 
quired intensity measurements (Fig[2]jc)). 

For each MZI, we characterised the phase by consider- 
ing the interferometer as an effective beam splitter with a 
reflectivity determined by the phase shift and the reflec- 
tivities of the four relevant beam splitters. This was again 
characterised using the ratiometric technique. Adjusting 
the voltage across the thermo-optic element varies the 
effective reflectivity, as shown in Fig. [2^d) . Fitting this 
data then provides both an estimate of the zero-voltage 
phase of the interferometer, and an independent consis- 
tency test of the four beam-splitter reflectivities which 
define each interferometer. These checks agreed within 
the measurement uncertainty. 



Demonstrating genuinely quantum operation and 
circuit complexity 

Having characterised the individual circuit elements 
classically, we now investigate the operation of the cir- 
cuit using quantum input light. We generated three 
individual input photons from two spectrally factorable 
down-conversion pair sources, using both photons from 
one pair and using the other as a heralded single pho- 
ton (see Methods for details). Detecting all four photons 
in this way to give four-fold coincidence measurements 
reduces the effects of noise terms. 

We first study the quantum interference from two- 
photon inputs to confirm our classical device charac- 
terisation and our single-photon indistinguishability. In 
these experiments, we injected photons into two of the 
selected input modes: cf, cd, or df (corresponding to 
regions I, II and III in Fig. |3ja)). We then measured 
interference visibilities for all possible coincidence out- 
comes by varying the timing of one photon using an off- 
chip optical delay stage (e.g. Fig. |3jd)-(e)), and com- 
pared the observed values to the quantum (boxes) and 
classical (triangles) predictions (see Fig. |3](a)). The ex- 
pected quantum visibilities were calculated directly from 
the quantum output state predicted by the classically 
characterised circuit unitary adjusted for the indepen- 
dently measured fidelity between photon pairs. The cor- 
responding classical visibilities were calculated by replac- 
ing the input Fock states with phase-averaged coherent 
states (see Methods). 

A chi-squared test verifies that these data are consis- 
tent with the quantum predictions whilst in strong dis- 
agreement with classical theory [IS]. The likelihood for 
observing a set of interference visibilities is calculated 
given measurement uncertainties in the observed inter- 
ference visibilities and the underlying circuit parameters 
(see Methods). The residuals from quantum theory give 
a reduced chi-squared Xr — 0.9; a value at least this 
large will occur with probability P(Xr > 0-9) = 0-5. By 
contrast, classical theory gives \ 2 r — 23 corresponding to 
P(Xr — 23) < 10~ 16 . Furthermore, the ultimate classi- 
cal visibility limit of 1/2 is expected and observed to be 
exceeded by output mode combinations de, cd, and dg. 

The three experiments each used different pair-wise 
combination of the three chosen input photons, allow- 
ing the different input photons to be tested individually 
against each of the other two. In all three cases, the 
experiments are in good agreement with the expected 
quantum results. This suggests that our element-wise 
characterisation technique is a good predictor of device 
performance and shows that each of the three photon 
quantum states has a good fidelity with the other two, a 
critical factor for observing genuine three-photon inter- 
ference. 

To demonstrate the complexity of the quantum cir- 
cuit, we study the higher-order nonclassical interference 
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FIG. 3. Two-photon interference (a) Experimentally measured two-photon interference visibilities (red circles) are 
compared against the quantum (clear boxes) and classical (blue triangles) predictions. Regions I, II and III contain experiments 
using input modes cf, cd, and df, respectively. The errors shown on the simulated quantum and classical visibilities were 
calculated by Monte-Carlo simulation; errors on the classical visibilities are smaller than the marker size, (b-c) The residuals 
between the measurements and the calculated classical and quantum visibilities, (d-e) Example two-fold coincidence counts 
between output channels dg and fg when two photons are input into modes df and the optical delay of mode d is varied. 



which arises when the three coupled interferometers are 
all operating simultaneously and in parallel, each injected 
with quantum light at the input. We observe this via the 
(heralded) three-photon coincidence counts, with three 
individual input photons coupled into modes c, d and 
/. Again detecting the fourth photon enabled discrimi- 
nation between downconversion events with one photon 
in each spatial mode (|lll) c( jf) from equally likely un- 
wanted noise events with two photons in each of only 
two input modes (|022) cdf ). After setting the temporal 
delay between input modes c and / to maximize their 
two-photon interference, we varied the delay for input 
mode d while simultaneously monitoring the eight three- 
fold coincidence combinations described in Fig. |4|a). We 
observed average four-photon coincidence rates of around 
16 mHz and measured the heralded three-folds continu- 
ously for 294 hours, iterating a full scan of the temporal 
delays each minute. This method averages out long-term 
systematic effects, such as drifts in the chip coupling ef- 
ficiencies and photon source performance [20], and al- 
lows an accurate calculation of the statistical error in the 
counts. 

As in the two-photon experiments, the observed three- 
photon quantum interference visibilities agree with quan- 
tum predictions based on the individually characterised 
circuit elements and are completely inconsistent with the 
equivalent classical predictions. The quantum prediction 
gives xl — 1-5 with P(xf > 1-5) = 0.2, whereas the clas- 
sical prediction gives Xr — 6 with P(Xr > 6) < 10~ 8 . 
Moreover, using a global optimisation routine, we deter- 
mined the maximum classical interference visibility for 
any possible circuit parameters with this circuit topol- 
ogy to be 0.59. This ultimate limit is exceeded by more 
than one standard deviation by output channel combi- 
nation cdf. Furthermore, a chi-squared test shows that 



the circuit parameters which result from this optimisa- 
tion are strongly inconsistent with our measured values 
(P(Xr > 4.4) = 10~ 8 ). Thus, only the full quantum ex- 
planation can plausibly account for the higher observed 
visibilities. We note that Fig. [4^a,) includes all mea- 
sured three-fold coincidence combinations, including sev- 
eral that occur very rarely and which therefore lead to 
large error bars in the measured visibility. Nevertheless, 
including all observed data in our analysis, we can ex- 
clude classical models with extremely high confidence, 
despite the low count rates in some channels. 

Finally, to verify that the observed interference re- 
sults from a quantum interaction of all three photons 
and is not explainable via a bi-separable interaction 
only, we also simulated the expected quantum visibilities 
when one of the three photons remains completely dis- 
tinguishable from the other two at all temporal delays. 
These quantum bi-separable explanations are also incon- 
sistent with the data, predicting at most a likelihood of 
P < 10~ 8 for the observed interference signature. 



DISCUSSION 

These results demonstrate genuinely multipartite 
quantum operation in a next-generation integrated cir- 
cuit which provides a critical new level of complexity in 
integrated quantum photonics. The measurements simul- 
taneously accessed three coupled interferometers, with 
classical interference at three circuit nodes and nonclas- 
sical interference at five circuit nodes. This experiment 
also represents the first observation of chip-based, mul- 
tipartite nonclassical interference which relies on more 
than two individual photonic inputs. The three-photon 
interferences were measured for eight different output 
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FIG. 4. Three-photon interference (a) The experimentally recorded three-photon interference visibilities (red circles) are 
compared against the predicted quantum (clear boxes) and classical (blue triangles) results as in Figure [3] (b-c) The residuals 
between our measured and simulated visibilities, (d-e) Example four-fold coincidence counts of output channels cdf and bdg 
when the input photon in mode d is temporally delayed. The shaded area shows the uncertainty in the determined visibility. 



combinations and are irreconcilable with classical mod- 
els. 

We develop a practical method to verify successful op- 
eration of the device: a loss-indcpcndcnt technique to 
classically characterise individual circuit parameters and 
two-photon interference to verify device performance un- 
derstood by simulations. Using a parameterised charac- 
terisation allows identification of poorly fabricated com- 
ponents and freedom to separately simulate individ- 
ual circuit subsections, such as on-chip state prepara- 
tion, manipulation and measurement. Having techniques 
which successfully predict device performance will be 
critical as experimental capacities continue to improve, 
making this demonstration an important step forward. 
An alternative scheme has recently been introduced for 
inferring the overall unitary transformation implemented 
by a device from a series of classical interference experi- 
ments using only the nominal input and output ports |21| - 
[23] . Whilst this does not give access to individual compo- 
nents, it may be useful for cases where only the operation 
of the device as a whole is of concern. 

In this paper, we highlight the often-unacknowledged 
fact that minimising losses will be critical for scaling up 
integrated circuits to the regime where they can no longer 
be simulated using classical processors. An array of on- 
going work seeks to do so by integrating [2"1H2"7] and syn- 
chronising [28 quantum light sources, as well as develop- 
ing high-efficiency integrated detectors [3D] ■ 

Increasing the complexity of integrated photonic de- 
vices requires not only an increased number of discrete 
optical modes but also complex multiphoton quantum 
interference across all of these modes. This work has al- 
ready verified the multiphoton interference necessary for 
the first nontrivial tests of recently proposed boson sam- 
pling problems [5TH55] . It also provides the first demon- 
stration of quantum operation of a chip which is suffi- 
ciently complex to allow a range of advanced quantum 



information protocols, such as teleportation and cluster- 
state generation, to be realised on an integrated platform. 



METHODS 

Device fabrication 

The waveguide circuit used in this work was fabricated by use of 
UV-direct write technology on a silica-on-silicon substrate [2]. The 
individual waveguides were written by focusing a continuous-wave 
UV laser (244 nm wavelength) onto the chip which is subsequently 
moved transversely to the surface normal with computer-interfaced 
2D motion control. The UV-writing process enables creation of 
beam splitters (X couplers) by crossing waveguides at different 
angles |34| . These are much smaller than traditional directional 
couplers which helps to reduce the effect of propagation loss in 
more complex circuits. The effective beam splitter reflectivity of 
these X couplers is primarily governed by the intersection angle 
of the guides which reduces sensitivity to wavelength, polarisation 
and temperature fluctuations making them extremely stable over 
the long experimental durations in this work. The thermo-optic 
phase shifters utilise a small NiCr electrode (0.35 fira X 50 fim X 
2.5 mm, 0.85 kOhm electrical resistance) deposited directly over one 
of the waveguides through which a current can be passed. The 
temperature-stabilised passive stability of the interferometers with 
the phase-shifters set to a constant voltage was measured to be less 
than 1° over 24 hours. 



High-brightness multiphoton states on-chip 

An 80 MHz Ti:Sapphire oscillator (Mai-Tai, Spectra Physics) 
produces 100 fs pulses at 830 nm (2.6 W average power) which are 
upconverted to 700 mW of 415 nm light in a 700 fim f} — BaB2C>4 
(BBO) crystal cut for type-I second-harmonic generation. This 
blue light is split on a 50:50 beam splitter and used to pump two 
8mm-long AR-coated Potassium Dihydrogen Phosphate (KDP) 
crystals phase-matched for degenerate type-II collinear paramet- 
ric down-conversion. We optimise the collection optics and spatial 
mode-matching to achieve a coincidence count rate of 160 kHz on 
each crystal with a raw heralding efficiency of 28-30 % without any 
filters. The source is designed to be spectrally factorable |35l which 
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improves the heralding efficiency we can achieve when interference 
filters (Scmrock, AA = 3 nm) are used to match the bandwidths 
of the broad and narrowband daughter photons. With the filters 
in place we achieve a four photon coincidence rate of 20 Hz and 
two-photon fidelities of 0.99 (narrowband-narrowband) and 0.96 
(narrowband-broadband) . 

Three of the photons were coupled into polarisation maintaining 
(PM) fibers and launched into the waveguide circuit using a butt- 
coupled PM v-groove array, fndex matching oil is applied between 
the fiber array and the chip to reduce reflection losses and a 6- 
axis piezo-controlled alignment stage provides all the degrees of 
freedom necessary to achieve optimal simultaneous coupling into 
all six input modes. The piezo-driven axes were operated in closed- 
loop mode to maintain this coupling throughout the experiment. 
An identical set-up is used on the output to a achieve maximal 
coupling from the chip into the single photon counting modules. A 
home-built FPGA-bascd logic unit records all desired coincidence 
counts simultaneously. 



Predicted visibilities 

The simulated visibilities in our interference experiments were 
calculated by first simulating the complete quantum output state, 
I '/'out), using the characterised circuit unitary, U c i rc . The intensity 
cross-correlation functions at zero and infinite temporal delay were 
then used to find an interference visibility. 

When photons are launched into modes {a;} of a linear optical 
circuit the intensity cross-correlation between output modes {bi} 
at zero temporal delay is: 



bi 



(i) 



where the intensity operator on mode i is Ij = b|h>i. Suppose 
the photon in mode undergoes a temporal delay then the total 
output state is now a classical mixture: 

po U t = i^r^X^uT^I + K^'X^I ■ (2) 

The intensity cross-correlation function at infinite delay is then, 
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and the visibility of the interference pattern between the n^-order 
output coincidence counts is then given by 
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The corresponding classical visibility of this interference pattern is 
given by injecting three equal amplitude coherent states of mutually 
randomised phase into modes {a;}. This ensures we mimic inde- 
pendent sources of light which will have no first-order correlation 
as required to compare against Hong-Ou-Mandel-type quantum in- 
terference. Coherent states represent the classical state which have 
the highest interference visibility, ensuring the bound we calculate 
is an upper limit. The resulting output vector of complex ampli- 
tudes is: 
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where e ou t is the vector of time-independent electric fields in each 



of the output modes and similarly for e; n 
intensity cross-correlation function is then 
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The classical cross-correlation function at infinite delay is calcu- 
lated in a similar manner, taking incoherent sums between the de- 
layed and non-delayed photons, 

V t} = (^ ■ £ E[ (\(4T* d) Kf + \(4lt\f) *>1- ■ ■ d0n 

(8) 

From which we calculate the classical interference visibility, 
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Chi-squared test 

To estimate the likelihood of observing a particular set of m in- 
terference visibilities v, we construct the probability density func- 
tion 
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where u>i(«) is the calculated i th visibility based on the set on n 
circuit parameters a, cr; is the measurement standard deviation in 
the observed visibility Vi, and Si is the measurement standard de- 
viation of the characterised circuit parameter ct^. We approximate 
each Wi(a) to be linear in a, 
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(11) 



which we verify to be accurate to within 0.02 over the range ±2sj. 
The resulting multidimensional Gaussian integral yields an analytic 
solution with an exponent quadratic in visibility residuals. We 
analyse our data in a basis Xj = Y_V Cij(vi — u>?) so that f(x) is 
factorable: 



/(x)~n cx p(- A ^?)' 



(12) 



The visibilities x arc thus statistically independent and a typical 
chi-squared test, \r = ^-/ m z\2T^i x i' can ^ e applied. If we in- 
stead errantly assume that the uncertainties in our simulated visi- 
bilities are uncorrelated, we calculate a chi-squared approximately 
0.3 lower than that reported for the quantum analysis of our ex- 
periments. 
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